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1. Molecule-based magnets 

Molecule-based magnets are systems whose principal structural building blocks are 
organic molecules lU |2]. We understand the microscopic origin of the magnetism 
that arises from such materials in terms of paramagnetic spin centres linked via 
superexchange pathways. The magnetic centres themselves may be atoms (such as 
s = 1/2 Cu^"*") or uncompensated spins on the molecules (as in the case of radical 
magnets) and the superexchange pathways will frequently pass through the molecular 
groups. The great promise of choosing to make our magnets out of molecules is 
tunability, that is, the possibility of manipulating the carbon chemistry of the molecules 
to form structures which promote magnetic exchange in certain geometries and to 
manipulate the strength of the coupling between the magnetic units. 

Research on molecular magnetic systems relies on a confluence of experimental 
and theoretical physics and chemistry. Chemistry is vital in the synthesis of new 
materials with the desired structural properties which incorporate the necessary 
magnetic ingredients; experiment allows us to probe the magnetism and we hope that 
theory can explain it. In this Comment we will discuss some examples of the use and 
successes of muon-spin relaxation (yU+SR) [3] as an experimental technique in this area. 
We will classify the systems that we discuss in terms of their magnetic dimensionality, 
that is, the number of dimensions in which the coupling between magnetic units is strong. 
The explanation of the physics of such systems is usually given in terms of low-energy 
models of exchange-coupled spins [H El E] whose long-wavelength properties have been 
elucidated using field-theoretic techniques and, in particular, the renormalization group 
(RG) O [8] . In presenting some possible future directions for muon spectroscopy in this 
field we will describe some of the theoretical models which are currently motivating the 
synthesis of new materials. 

The subject of molecular magnetism has been reviewed thoroughly elsewhere [HE], 
as has early, and more recent, work using /i+SR P |TOl [11] [121 113]. However, it is worth 
mentioning some of the notable themes from the early muon work, which include: (i) the 
observation of long-range magnetic order (LRO), such as in the purely organic material 
/9- phase ]9-NPNN, where oscillations in the muon asymmetry are found below 0.7 K |14j : 
(ii) the observation of dynamics in molecular nanomagnets such as the single molecule 
magnet Mni2Ac [15] and (iii) the possibility of observing collective magnetic behaviour 
such as the spin-crossover transition [9J. These investigations play to the strengths 
of /i"'"SR in that they require a probe sensitive to small magnetic moments and their 
dynamics and the possibility of distinguishing between uniform magnetic effects and 
minority, strongly magnetic phases. 

The Comment is structured as follows: in section [2] we will discuss the theoretical 
prediction and realization of one-dimensional magnetism from molecular chains. In 
section [3] we turn to layered two-dimensional systems and in section [H to the recently 
realised possibility of switching between dimensionalities. Finally we discuss some of 
the possible future directions for this field in section [51 
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Figure 1. The structure of Cu(pyz)(N03)2- (a) The Cu-pyz-Cu chains, (b) The chain 
packing viewed along the a-direction. Figure taken from Ref. [50] . 

2. Magnetism in one dimension 

Much recent work on molecular magnets using //"""SR has concerned materials whose 
low-energy physics may be thought of as arising from interactions which are constrained 
to act in less than three spatial dimensions. Such physics usually achieved by forming 
an array of transition metal ions, such as s = 1/2 Cu^"*" and linking these using one 
or more molecular ligands in such a way as to promote strong interactions between the 
magnetic ions only along particular directions and weaker ones along others. 

As an example, we will consider a material whose structure promotes strong, 
antiferromagnetic superexchange along one-dimensional, parallel, linear chains, and 
weaker exchange in the other two perpendicular directions. Such a material may be 
constructed if we link Cu^+ ions along one direction using pyrazine (pyz = C4H4N2) 
ligands and separate the chains using NO3 groups. The resulting material formed from 
this recipe is Cu(pyz)(N03)2 and is shown in figure [H 

To understand the physics of such materials we seek to describe them using simple 
spin models. To a good first approximation the crystal field surrounding the Cu^^ ions 
does not cause the spin to favour any particular spatial direction, resulting in a so- 
called Heisenberg spin [16]. The resulting one-dimensional {d = 1) quantum Heisenberg 
antiferromagnet (IDQHAF) may be described by a model Hamiltonian [U |5l [17] 

H = JJ2S.-S„ (1) 

(ij) 
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Figure 2. (i) The spinon excitation in a d = 1 chain, (ii) The phase diagram of a d = 1 
chain system with nearest neighbour couphng Ji, next- nearest neighbour coupling 
J2 and anisotropy parameter a [VBS=valence bond solid; T-L liquid=Tomonaga 
Luttinger liquid (see Section [5|).] (iii) The valence bond solid with a single s = 1/2 
excitation, (iv) The phase diagram of a d = 2 system showing a quantum critical point 
at a coupling and a quantum critical region at elevated temperatures. 



where i and j label spins along the same chain and the sum is over unique nearest 
neighbours bonds only. For materials such as Cu(pyz)(N03)2 the parameter J may be 
determined from the result of magnetic susceptibility or pulsed field magnetization (see 
below) . 

The model Hamiltonian in eqn[T]is very interesting and has motivated a wide range 
of research [HIS]. It possesses a continuous 0(3) symmetry (here the global rotation of 
all of the spins through an arbitrary angle). In common with other d = 1 models with 
finite range interactions there will be no LRO for T > 0. (For d = 1 this is easy to see: 
pass the message "line-up" along the chain in the presence of a defect (or a fluctuation) 
and the message is lost at the weak link: since each spin is only able to communicate with 
its neighbours the message cannot pass around the weak link). Although this problem 
may be solved exactly using the famous Bethe ansatz [6] the solution is rather unwieldy 
and many of our insights arise from the predictions of a continuum model which, as 
long as the long-range spin-spin coupling is not too large (see section [5]), results from 
mapping the problem onto a d = 1 model of spinless fermions and then using a technique 
known as bosonization. The predicted state is known as a Tomonaga-Luttinger (T-L) 
liquid [5l [6] and its excitation spectrum is one of gapless, linearly dispersing, particle- 
like excitations. Some physical insight into the nature of these excitations comes from 
figure El^i). An antiferromagnetically correlated region [as shown in (a)] may be excited 
if we flip a section of the spins, as has been done after the fourth spin in (b). This 
creates a spin 1/2 excitation, known as a spinon [shown shaded in figure [2]^i)], which 
may move by hopping by two lattice spacings at a time [shown in (c)]. 
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Figure 3. Zero field muon-spin relaxation data measured on Cu(pyz)(N03)2, showing 
long-range magnetic order below Tn = 107 mK [21]. . 



The idealisation that leads to the T-L model has ignored many interactions in 
Cu(pyz)(N03)2 and one of these proves particularly relevant: the interaction between 
chains parametrised by an exchange interaction J^. This changes the Hamiltonian to 

where i and j' label spins on adjacent chains. In order to probe d = 1 behavior a 
system described by this model we would need to work in a temperature regime where 
<^T <^ J . Once we start to approach a regime where the characteristic energy of the 
thermal fluctuations allows detail on the level of to be resolved, we should expect that 
the three-dimensional nature of the material becomes evident and it will magnetically 
order. However, the strongly anisotropic nature of the exchange {J±/ J -C 1) will have 
consequences on the ordered state that arises [18] and the ordering temperature Tn is 
driven downwards by the one-dimensional nature of the system compared to a three- 
dimensional one. The ordered magnetic moments will also be renormalized, making 
them potentially very small, and difficult to observe with magnetic susceptibility. Finite, 
but possibly very long, correlation lengths may exist above the ordering temperature, 
reducing the entropy in the system. Consequently, the entropy change on ordering 
will be much reduced compared to a more isotropic system and may prevent specific 
heat measurements from detecting a transition. Add the difficulty of deuterating all 
of the hydrogen- containing groups in the system (and the possibility of deuteration 
changing the delicate superexchange pathways [I9]) and neutron diffiaction also becomes 
difficult. It is here that /i+SR has proven very effective in detecting magnetic transitions 
in reduced dimensions that have proven difficult to observe using more conventional 
experimental methods. Before the /i^SR study, no magnetic order had been observed 
in Cu(pyz)(N03)2 despite specific heat measurements made down to 70 mK. However, 
the clear zero-field precession signal from /i+SR [21] shows unambiguous evidence for 
LRO below a transition temperature of Tn = 0.107 K as shown in figure [31 

A useful figure of merit in comparing the success with which a material realizes 
one- dimensionality is Tn/J, since this quantity should be zero in the ideal case and 
of order and close to unity for an isotropic material. This quantity may also be used 
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to estimate J± , which is the parameter of interest in the Hamiltonian in eqn [2J One 
method of doing this is has been developed from Quantum Monte Carlo computations 
[22], whose results may be approximated using the expression 




with a = 2.6 and c = 0.233. This yields a value of \J±/J\ = x lO^^ for 
Cu(pyz)(N03)2, demonstrating that it is indeed a very well isolated approximation of 
a (i = 1 magnet. In addition, coupled chain mean field theory [23] predicts an ordered 
moment of 0.05 /xb- 

The application of a magnetic field to our d = \ system acts to polarize the spins. 
If we imagine doing this at T = then there will be a point at a field [i^B ^ J where all 
of the spins will be aligned. This aligned state is separated from the partially polarised 
antiferromagnet by a quantum critical point (QCP) [5] at a field = 2J/g^B (we 
neglect interchain interactions here). One advantage of investigating molecular magnetic 
systems such as these is that the exchange strength, and hence the QCP, is typically of 
order J fa 10 K, translating as ~ 10 T. Although such fields are currently outside 
the range accessible via /i"'"SR they are accessible using pulsed magnetic fields and the 
combination of high- field magnetization and muon spectroscopy [23] has proved valuable 
in fully characterising new materials in recent years. These characteristic energy scales 
are orders of magnitude smaller than the corresponding ones for inorganic spin chains 
such as Sr2Cu03, whose QCP occurs at an experimentally inaccessible field around 200 
times larger than that of Cu(pyz)(N03)2. 

An even more ideal example of this physics is found in the radical system DEOCC- 
TCNQF4 [25]. In this case no LRO transition is observed using /i"'"SR down to 20 mK, 
constraining \J±/J\ < 6 x 10~^, making this the most successful realization of a 
IDQHAFM yet reported. In this case, the weak relaxation that is observed via yU"'"SR 
allows the investigation of the magnetic fiuctuation spectrum. Whether the excitations 
in IDQHAFs are ballistic or diffusive remains a key question in the physics of spin chains 
[26l [27] and one that can be addressed with LF /i+SR measurements. Field-dependent 
studies can be used to distinguish between these types of spin transport, since their 
spin autocorrelation functions have different associated spectral densities f{uj) which 
vary as for d = 1 diffusive transport and f{ui) oc ln(J/w) for ballistic motion. In 
DEOCC-TCNQF4 diffusive transport was identified. This contrasts with the inorganic 
IDQHAF Rb4Cu(Mo04)3 where ballistic transport was shown to dominate the muon 
relaxation [21]. The details and general rules of the d = 1 transport are complicated 
and this continues to be an area exciting much interest and new research. 

3. Two dimensions 





(3) 




The model of the two-dimensional square lattice quantum Heisenberg antiferromagnet 
(2DSLQIIA) has received significant attention as another apparently simple system 
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Figure 4. Zero field ^+SR data measured on [CuHF2(pyz)2](BF4). (a) Zero field 
oscillations below Tn = 1.54 K allow us to identify LRO and probe the critical 
behaviour. [S^ (b) For T > Tn we observe muon-fluorine dipole-dipole oscillations 
that allow us to deduce muon sites [33l [34] . 



where the influence of dimensionality leads to the emergence of an unusual ground state 
and excitation spectrum [29]. It is also, of course, the low-energy model that describes 
the Mott insulating parent phase of the cuprate superconductors. The 2DSLQHA shows 
no magnetic order for T > 0, owing to the Coleman- Mermin- Wagner theorem [30l 131] . 
but there is good evidence that it orders at T = with a much renormalized ordered 
moment of around 0.62 /xb [29]. The excitations of the model at T = are gapless spin 
wave excitations with a linear dispersion [H El |29] • 

Of the materials recently synthesised which show two-dimensional {d = 2) magnetic 
properties, one of the most successful classes is based on two-dimensional square layers 
of Cu^"*" linked with pyrazine groups to form [Cu(pyz)2]^^ layers with an intralayer 
exchange coupling J. These square layers may be linked with HF^~ ions to form 
pseudocubic arrays of [CuHF2(pyz)2]X [shown inset in Fig. Hj^a) [32]], where X is an 
anion that sits in the centre of the cube and provides charge neutrality. As in the d = 1 
case, we are able to identify magnetic order in this system at low temperature and assess 
the degree to which it may be described as two-dimensional [33] . As an example we may 
consider [CuHF2(pyz)2]BF4 (Fig. H]) which magnetically orders at T = 1.54 K [32]. If we 
attribute the ordering to an interlayer coupling J± we may use the results of quantum 
Monte Carlo calculations which, in this case, predict [22] that 

l^^/^l=exp(6-^,), (4) 
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Figure 5. (a,b) Transverse field /i+SR measurements on Cu(HF2)(pyz)2]BF4. (c) The 
relaxation rate allow us to extract the ordering temperature in field, (d) The result is 
a non-monotonic phase boundary caused by a small XK-like component in the spin 
Hamiltonian of the system '33' . 



where ps = 0.183 J is the spin stiffness and the parameter b = 2.43. For 
[CuHF2(pyz)2]BF4 we find J = 6.3 K and \ J±/J\ = 9 x 10~^, making this material 
a highly anisotropic d = 2 antiferro magnet. 

In addition, the fact that these systems contain fluorine leads to oscillations at 
temperatures above Tn attributable to F-yU"*" dipole-dipole coupling [Figure Hl^b)]. These 
are useful here as they allow us to work out where the muons giving rise to these signals 
are localised with respect to the fluorine nuclei in these materials [M]- It is worth noting 
however that in complex materials such as these there are often several classes of muon 
sites, not all of which may be accounted for by those muons giving F-yU"*" oscillations. 

Although interlayer coupling in a 2DQIIAF again provides a means for magnetic 
order to occur at non-zero temperatures there is evidence that in this system a degree 
of ^f-factor anisotropy causes the Cu^"*" ions to take on a degree of XF-like behavior 
[35] . This is of interest as the pure XF- model allows the occurrence of a new sort 
of excitation: the vortex. A single vortex is unstable alone, and in the absence of 
an external electromagnetic field, we may only stabilise vortices at low temperature 
by binding them with antivortices. The transition between bound vortex-antivortex 
pairs and single vortices stabilised by thermal fluctuations at high temperature was 
originally identified by Kosterlitz and Thouless and by Berezinskii for an ideal system, 
but has never been observed in a bulk magnet [Ml [371 EH] . However, the physics that 
those authors identified may have an effect even when only a small XF-like component 
is present. In this case we expect that, while at high temperature the spins will be 
Heisenberg-like, as the temperature is lowered to a scale where the XF-like component 
is resolvable it becomes energetically favourable for the spins to tip into the XF-plane. 
On further decrease of temperature the possibility of a vortex-binding transition greatly 
increases the correlation length, at which point the interlayer coupling J± causes the 
system to undergo a transition to LRO [35] . 

One experimental consequence of XF-like behaviour is that the phase boundary 
in applied magnetic field takes on a non-monotonic form as shown in figure [5] [39l |33] . 
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The application of the field makes it energetically favourable for the spins to tip into 
the XY plane, effectively reducing the dimensionality. It then has a dual effect on 
the spins by suppressing the amplitude of the order parameter by forcing spins along 
a particular direction and also reducing the phase space for phase fiuctuations from a 
sphere to a circle. At low fields the magnetic ordering transition is controlled by phase 
fiuctuations and at higher field by amplitude fiuctuations leading to the non-monotonic 
phase boundary that is observed in experiment [39] . We have been able to measure the 
phase boundary in field using both LF /x+SR and transverse field (TF) /i+SR [33] as 
shown for [Cu(HF2)(pyz)2]BF4 in figured It has been argued [39j that the energy scales 
for this physics are controlled by the Kosterlitz-Thouless physics described above. It is 
probable that the physics of vortices, along more complicated topological objects such 
as merons and skyrmions, will be an area of future interest in this field. 

4. Switching the dimensionality 

We might ask if it is possible to gain control over the dimensionality of a system. There 
have been two recent advances that suggest that this is a fruitful area for future work: the 
first a synthetic route; the second using pressure as a means of creating a dimensionality 
switch. 

The first route [40j relies on varying the proportions of reactants in the 
material synthesis, which makes it possible to create two distinct, but closely related 
materials, one two-dimensional; the other one-dimensional. The coordination polymer 
[Cu(pyz)2(pyO)2](PF6)2 (where pyO denotes pyridine- N-oxide=C5H5NO) is another 
d = 2 square lattice material based on [Cu(pyz)2]^"'" layers, this time with pyO groups 
protruding from them. By changing the proportion of pyz and pyO molecules in the 
preparation from 3:1 to 2:1 the water- containing material [Cu(pyz)(pyO)(H20)2](PF6)2 
is produced. Here pyz ligands link Cu atoms along the 6-axis only, forming highly 
one-dimensional chains. Crucially, the strong nearest-neighbour exchange interaction 
maintains a value of J ~ 8 K in both cases. 

These structures suggest d = 2 and d = 1 magnetic behaviour, which is indeed 
confirmed by a combination of high- field magnetometry and muon-spin relaxation [ID] . 
The results of the /x"^SR measurements are shown in figure [6] where we are able to 
establish where each material undergoes a transition to LRO. The ordering temperature 
and magnetic moment size of the d = 2 material are both larger than in the analogous 
d = 1 material, which refiects the expected increase in the infiuence of fiuctuations for 
d = 1 compared to = 2. We conclude, in this case, that to turn a material from d = 2 
to d = 1, one must just add water! 

The second method of controlling dimensionality involves using pressure as a switch 
[20] . The material CuF2(H20)2pyz [Hj is based on a chain-like structure but, rather 
counter-intuitively, it is a good realization of a 2DQHAFM. The reason for this is that 
the chains arrange themselves in such a way as to form square layers of Cu^^ ions with 
superexchange pathways provided by hydrogen bonds. Crucially the unfilled d^^-yi 




orbital lies in the plane of the square layers, providing strong two-dimensional exchange 
with only weak coupling along the chains. 

Applying hydrostatic pressure causes a dramatic change in behaviour above around 
10 kbar [20]. Here the magnetization shifts from that typical of a d = 2 system to that 
of a c? = 1 magnet. /i"'"SR has been carried out under pressure and reveals the change 
in Tn as expected from a change from d = 2 to d = 1. The microscopic explanation for 
this quantum phase transition is suggested by a combination of x-ray diffraction and 
electron spin resonance [l2lll3]. The application of large hydrostatic pressures causes a 
small structural distortion which, in turn, causes the Jahn- Teller axis of the Cu^"*" ions 
to switch directions, resulting in the d^2_y2 orbital shifting from linking Cu atoms in 
hydrogen bonded [Cu(pyz)2] layers to point along the Cu-pyz-Cu chains. 

It is possible that such a switching of magnetic dimensionality will be found, and 
ultimately engineered, in a variety of magnetic coordination polymers. The fundamental 
change in the character of a system in going from ci = 2 to c? = 1 is such that this 
switching feature could conceivably find a use in future quantum technologies. 

5. Frontiers: a zoology of ground states 

We now turn to a possible future direction for the engineering of coordination polymers. 
We hope that it will soon be possible to produce materials that realise still more 
complicated models of low- dimensional magnetic behaviour. These would allow us to 
realise a large number of new and potentially exotic ground states. Here we briefly 
review the predictions for some of the quantum states that may be achieved if we are 
able to tune the interaction strength of our s = 1/2 spins. For the full story we refer 
you to some of the excellent books on this subject [U El E]. 

5.1. d=l: Tomonag a- hutting er, Valence bond solid and Ising-N eel states 

When describing the IDQHAFM in section[2]we claimed that the disordered T-L ground 
state with its gapless linear excitation spectrum was reached if the coupling was not 
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too large. More precisely, we mean that that the antiferromagnetic nearest neighbour 
exchange coupling Ji is dominant and the antiferromagnetic next-nearest neighbour 
exchange J2 is small in comparison. However, if the ratio of the couplings exceeds a 
critical value J2/J2 > Jc then the system adopts another ground state: the spin-Peierls 
or (in more modern terminology) the valence bond solid (VBS) state, illustrated in 
figure [2|^iii). The VBS ground state involves the spins staggering their bond exchange 
energy, essentially forming antiferromagnetically coupled dimers along the chain. There 
are two ways of doing this along any chain (since a particular site may form a dimer with 
its neighbour to the left, or to the right) and so the ground state that is realised must 
break translational symmetry. An antiferromagnetically coupled dimer has a singlet 
ground state, resulting in a total S* = ground state for the VBS. The excitations of 
the VBS phase are lone spins sat between the valence bonds. These are spin s = 1/2 
excitations which, as in the T-L case are also known as spinous but here there is an 
energy gap A which must be exceeded to create such an excitation. In fact, the energy 
required to excite a spinon of momentum p is given by Ep = (A^ + fpP^) ^ , where vp is 
the Fermi energy of the system. 

Although we have hitherto concentrated only on low-energy models which include 
purely spin degrees of freedom, the realization of the VBS state also requires the 
inclusion of magnetoelastic coupling (and it is in this context that it is usually called 
a spin-Peierls state). Above the spin-Peierls transition temperature Tsp, we have a 
uniform antiferromagnetic next-neighbour exchange in each chain; below T^p there 
is an elastic distortion resulting in the dimerization described abve, and hence two, 
unequal, alternating exchange constants. The alternating chain possesses an energy 
gap between the singlet ground state and the lowest lying band of triplet excited 
states which closes up above Tsp. The transition temperatures may be related to the 
relevant coupling constants; whereas the conventional Peierls distortion is expected at 
a temperature Tp ~ (i^p/^B) exp(— 1/A), where Ep is the Fermi energy of the system 
and A is the electron-phonon coupling constant, the spin-Peierls transition is expected 
at TsP ~ (J/fcp) exp(-l/A). J <. Ep and hence T^p -C Tp [B]. 

One example system is the organic salt MEM(TCNQ)2, which consists of one- 
dimensional stacks of planar TCNQ molecules, each of which has a charge of —\e 
associated with it. Adjacent stacks are separated by arrangements of MEM molecules, 
each of which possess a localized charge of +e. It undergoes two structural distortions. 
The first, which occurs at 335 K, is a conventional Peierls transition in which the TCNQ 
chains dimerize. This results in a change from metallic to insulating behaviour as a 
single electronic charge becomes localized on each TCNQ dimer; the single spin on each 
dimer couples antiferromagnetically to its neighbours. This phase persists down to the 
spin-Peierls transition at 18 K, where a further dimerization of the TCNQ stacks takes 
place (this is a tetramerization of the original chain). /iSR studies indicate a slowing 
down of the electronic fluctuations resulting from the opening of the spin-Peierls gap 
in the magnetic excitation spectrum as the temperature is lowered below T^p [15]. At 
the very lowest temperatures the electronic spin fluctuations freeze out and the muon- 
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spin depolarization is dominated by persistent slow fluctuations which can be ascribed 
to a defect-spin system [16]. It has been speculated that the muon spin might locally 
create a spin defect by breaking a singlet pair, giving rise to the relaxation ascribed to 
defect spins and which is only revealed when the other sources of relaxation have frozen 
out [IE]. These results have now been extended to the related material DEM(TCNQ)2 
|47] . Other examples of spin-Peierls systems include other organic systems such as 
TTF-CuS4C4(CF3)4 (T,p = 12 K) [48j and TTF-AuS4C4(CF3)4 {T,p = 2 K) [Ml SI]. 
A common feature of such materials is flat organic molecules in columnar stacks. The 
large interchain separation and weak van der Waals intermolecular interactions favour 
the dominance of magnetoelastic effects over interchain ordering. In contrast the chains 
in corresponding inorganic materials, such as copper chain compounds, are quite rigid 
due to the ionic bonding and only a single example of an inorganic spin-Peierls material 
is known (CuGeOa, with T,p = 14 K ^SOJ). 

As the previous example of the VBS illustrates, allowing strong interactions from 
distant spins opens up new possibihties of interesting phases to study. In fact a rich seam 
of physics is revealed if start with a. d = 1 spin Hamiltonian, and allow the possibility of 
both an XF-like spin anisotropy component and next-nearest neighbour interactions. 
This is the content of the J1-J2 model, whose Hamiltonian is given by 

H = J, Y^iStS] + SfS] + aStS^) + J, Y.{StSl + SfSl + aStSl), (5) 

{ij) ik 

where i and j label nearest-neighbours; i and k label next-nearest neighbours and a 
is a parameter which we may vary from {XY spins), through 1 (Heisenberg spins) 
to large values (resulting in Ising spins). This model may, in turn, be mapped onto 
the sine-Gordon model of quantum fleld theory [5l [51] (the origin of the latter's curious 
name is described in Ref [51] ) and whose RG scaling analysis leads to the phase diagram 
shown in flgure [2]^ii). 

From flgure|2t^ii) we see that tuning the parameters of this model allows us to access 
three ground states: the T-L liquid; the VBS (both described above) and the Ising-Neel 
state. [It is worth noting that along the line a = 1 (corresponding to Heisenberg spins) 
it can be argued that for J2/J1 < Jc the system is a T-L liquid but for J2/ Ji > Jc 
the system adopts the VBS state (as claimed above and in section [2])]. The Ising-Neel 
state, reached by making a > 1, with moderate next-nearest neighbour coupling, is 
characterized by a spontaneously broken symmetry at T = 0. Here we have magnetic 
LRO (not seen in either the T-L or VBS states) involving a staggered expectation value 
for the z-component of the spins. The lowest energy excitations of this phase are gapped 
s = 1/2 spinous. The properties of these states are summarised in table [TJ 

5.2. d = 2: quantum criticality and spin liquids 

The ground state for the 2DQHAF is antiferromagnetically (Neel) ordered at T = 
[U [5l |29]. This phase has gapless spin- wave excitations which do not carry a deflnite 
spin. (However, if the Neel order is aligned along z, the spin waves have total Sz = ±1.) 
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Phase 


T = ground state 


excitation 


character 


Tomonaga-Luttinger liquid 


Spin disordered 


s = 1/2 spinous 


gapless 


Valence bond solid 


Dimerised 


spinous 


gapped 


Ising-Neel 


Neel ordered 


spinous 


gapped 



Table 1. Phases of the d ~ 1 magnet. 



Phase 


T = ground state 


excitation 


character 


Neel {g < g^) 


Staggered order 


Total Sz = 1 spin waves 


gapless 


Quantum paramagnet {g > gc) 


disordered 


s = 1 flipped spins 


gapped 


Valence bond solid {g > gc) 


total S = 


s = 1 


gapped 



Table 2. Phases of the d — 2 Heisenberg magnet. The inclusion of Berry phases in 
the theoretical description leads to the VBS phase for g > gc- 



As we increase the strength of the effective long-range coupling (which we call g here) 
another phase is realised above a QCP at a critical effective coupling gc, as shown in 
figure El^iv). This phase, which is often called a quantum paramagnet, does not show 
LRO and has a gap against s = 1 excitations, which may be thought of as flipped 
spins. This state can also be realised in the sort of spin model we have been considering 
by including next-nearest neighbour interactions as we discussed for the d = 1 model. 
In fact, numerical and series expansion studies suggest that gc corresponds to a ratio 
between next-nearest and nearest neighbours of J2/J1 = Jc ~ 0.4. The arguments 
leading to the quantum paramagnet neglect the contribution of Berry phases to the 
final configuration of the system. If these are included a VBS state is realised with 
singlet bonds lining up in columns or plaquettes on the d = 2 square lattice. Such a 
ground state is four-fold degenerate and so, as in the d = 1 case, the VBS state that 
is realised must break translational symmetry. However, the lowest lying excitation is 
the same as predicted for the quantum paramagnet (a gapped s = 1 quasiparticle), but 
higher excitations will be different [5]. 

Finally, although everything we have discussed in this section applied for T = 0, 
we may consider non-zero temperatures in this d = 2 case, which results in the phase 
diagram shown in figure [2]^ iv). The curved lines represent crossovers in behavior which 
allow us to identify a quantum critical region, which extends down to T = for 
g = gc- Here, both quantum fluctuations (i.e. time dependence in the wavefunctions 
deriving from uncertainty relations) and thermal fluctuations are equally important 
and it is impossible to describe the excitations in terms of nearly-free quasiparticles. 
In this region the phase of the wavefunctions are as incoherent as Nature allows, 
realising an (almost) ideal fluid. Quantum critical behaviour continues to be an area 
of much experimental and theoretical interest. The properties of the d = 2 system are 
summarised in table [2] 

As an example of how muons may be useful in elucidating the physics of the exotic 
phases described here, we mention the case of k-(BEDT-TTF)2Cu2(CN)3 [52j. This 
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(a) (b) (c) 

Anisotropic triangular lattice ID (decoupled chains) 2D {'square' lattice) 




Figure 7. (a) The triangular lattice with couplings ti, t2 and ^3 as shown. By varying 
the coupling we may realize (b) the d = 1 or (c) the d — 2 cases, (d) Generalized phase 
diagram [S3] for the triangular lattuce. (QSL=quantum spin liquid phase, k-CN = k- 
(BEDT-TTF)2Cu2(CN)3.) (e) Phase diagram of k-(BEDT-TTF)2Cu2(CN)3 derived 
from /x"*"SR measurements |52j . 

is a c? = 2 system based around a triangular lattice of (antiferromagnetically coupled) 
s = 1/2 spins, whose low energy physics provides an example of a frustrated, non- 
collinear antiferromagnet. In fact, the triangular lattice, shown in figure Wi^) may be 
viewed as a meeting point of the d = 1 and d = 2 cases. We may recover those cases by 
tuning the coupling constants to form either well coupled d = 1 chains (e.g. by taking 
ti ^ ^2,^3) [figure El^b)] or, with the choice ti = ^2 ^ ^3, forming a lattice which is 
topologically equivalent to the d = 2 square lattice case [figure [Tl^c)]. A general phase 
diagram for the triangular lattice, derived from calculations using modified spin wave 
theory [53] and assuming perfectly localized Heisenberg spins, is shown in Fig. Ej^d). 
The d = 1 case is realized far along the horizontal axis, where the T-L ground state is 
identified as one example of a quantum spin liquid (QSL) state (here used as a general 
term used to describe a T = ground state without LRO, of which there are potentially 
very many varieties). The d = 2 case is found far along the vertical axis showing the 
expected T = collinear (Neel) LRO. Of note in figure [Tj^d) is the isotropically coupled 
case (ti = t2 = ts) where the classical solution of antiferromagnetic LRO with a 120° 
spin structure is realised, albeit with moments much reduced by quantum fluctuation. 
This neighbours a QSL state for the case where the couplings are slightly anisotropic. 
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Returning to the case of k-(BEDT-TTF)2Cu2(CN)3, we see that its couphngs 
appear to place it on the border of the QSL and Neel phases in figure[7|^d). In fact, a more 
reahstic description of the material is acheived by assuming that U/t (the ratio of the 
electronic correlation energy to the bandwidth) is finite. This, which amounts physically 
to a relaxation of the assumption that the relevant electronic degrees of freedom are 
perfect Heisenberg spins (or so-called "quantum tops" ) is expected to destabilise phases 
showing LRO near the origin in figure [Tl^d), enlarging the QSL phase to such an extent 
that K-(BEDT-TTF)2Cu2(CN)3 should be expected to realize a QSL ground state. In 
fact, the /i"'"SR measurements [52] suggest that the material realizes a Z2 QSL ground 
state showing topological order [5], with total S = and s = 1/2 spinon excitations 
with a very small gap against their creation. (The "Z2" assignment tells us that the 
system is described by a topological gauge field theory |51] somewhat similar to those 
used to describe fractional quantum Hall fluids, whose description is locally invariant 
under the Z2 group of transformations of the fields.) There is also evidence for another 
type of gapped excitation in this system, which is a point-like vortex known as a vison. 
For fi;-(BEDT-TTF)2Cu2(CN)3 muon measurements |52] not only suggest this complex 
excitation spectrum, but the application of magnetic fields allows the system to be 
driven into other phases showing LRO with extremely small magnetic moments, leading 
to the phase diagram shown in figure [Tl^e). The exquisite sensitivity of /i'^SR and the 
potential to apply it to samples which are smaller than those required for inelastic 
neutron measurements have made it particularly useful in this context. 

6. Conclusion 

In the short term muons will continue to allow us to probe low-moment LRO and aspects 
of the excitation spectrum in low- dimensional magnetic systems based on molecular 
building blocks. In addition to the cases outlined here there are many other avenues 
for future study in this field, including the investigation oi d = Q systems such as 
those based on dimers or on molecular nanomagnets formed from clusters containing 
large numbers of strongly coupled spins. Also of interest are spin-ladder systems which 
interpolate between the harsh extremes of ci = 1 and d = 2. In the future, the possibility 
of being able to work in more extreme conditions involving large applied magnetic fields 
and pressures will open out more of the phase diagrams of these systems. At the same 
time, advances in materials synthesis will allow us to approach the platonic ideals of the 
models described in this Comment. 
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